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ABSTRACT 

We present maps of the plane-of-sky magnetic field within two regions of the Taurus molecular 
cloud: one in the dense core L1495/B213 filament, the other in a diffuse region to the west. The 
field is measured from the polarization of background starlight seen through the cloud. In total, 
we measured 287 high-quality near-infrared polarization vectors in these regions. In L1495/B213, 
the percent polarization increases with column density up to Ay ~ 9 mag, the limits of our data. 
The Radiative Torques model for grain alignment can explain this behavior, but models that invoke 
turbulence are inconsistent with the data. We also combine our data with published optical and 
near-infrared polarization measurements in Taurus. Using this large sample, we estimate the strength 
of the plane-of-sky component of the magnetic field in nine subregions. This estimation is done with 
two different techniques that use the observed dispersion in polarization angles. Our values range 
from 5 — 82 /iG and tend to be higher in denser regions. In all subregions, the critical index of the 
mass-to-magnetic flux ratio is sub-unity, implying that Taurus is magnetically supported on large 
scales (~ 2 pc). Within the region observed, the B213 filament makes a sharp turn to the north and 
the direction of the magnetic field also takes a sharp turn, switching from being perpendicular to 
the filament to becoming parallel. This behavior can be understood if we are observing the rim of a 
bubble. We argue that it has resulted from a supernova remnant associated with a recently discovered 
nearby gamma-ray pulsar. 

Subject headings: dust, extinction — ISM: bubbles — ISM: clouds — ISM: individual (Taurus) — ISM: 
magnetic fields — Polarization 



1. INTRODUCTION 

It has been known for decades that observed starlight is 
polarized (Hall 1949; Hiltner 1949). This polarization is 
generally understood to be caused by non-spherical dust 
grains between Earth and the star that preferentially 
align with their long axes perpendicular to the local mag- 
netic field direction. At optical to near-infrared wave- 
lengths, the grains will preferentially absorb starlight 
that is linearly polarized parallel to their long axes. 
Thus, it is possible to probe the interstellar magnetic 
field by observing the polarized light because at optical 
and near-infrared wavelengths, the measured polariza- 
tion angle will be parallel to the magnetic field direction. 

The importance of magnetic fields to star formation is 
a hotly debated topic. Many to most molecular clouds 
are gravitationally bound, yet the observed star forma- 
tion efficiency of clouds is only a few percent (Myers 
et al. 1986). Magnetic fields are one mechanism that 
can provide cloud support against gravitational collapse 
(Mouschovias 1976). An alternative theory is that mag- 
netic fields are too weak to resist gravity and it is in- 
terstellar turbulence that regulates star formation (Mac 
Low & Klessen 2004). 

The Taurus molecular cloud is an excellent target for 
exploring the importance of magnetic fields in molecu- 
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lar clouds. As one of the closest low-mass star-forming 
clouds at a distance of 140pc (Elias 1978; Kenyon et al. 
1994; Wichmann et al. 1998; Loinard et al. 2005, 2007; 
Torres et al. 2007), Taurus has been a frequent target 
for studies (see the extensive review by Kenyon et al. 
2008). Furthermore, recent work has revealed a coupling 
between the gas and the magnetic field in Taurus. Gold- 
smith et al. (2008) found striations in their 12 CO data 
that matched the angle of the optical polarization vec- 
tors measured by Moneti et al. (1984). Goldsmith et 
al. offered two possible mechanisms to explain the stri- 
ations and their observed alignment with the magnetic 
field. Large-scale flows are seen in the 12 CO and 13 CO 
channel maps; given the coupled magnetic field and gas, 
shear in these flows could explain the striations. Alter- 
natively, a magnetosonic wave traveling perpendicular to 
the field could compress or rarefy the gas into the per- 
ceived striations. The filamentary region B213 in Taurus 
also has an apparent connection between the gas and the 
magnetic field. The optical and near-infrared polariza- 
tion vectors are oriented perpendicular to the long axis 
of the filament (Heyer et al. 1987; Goodman et al. 1992) 
implying that gravitational collapse has occurred along 
the field lines but not across them. 

Several studies have measured the large-scale magnetic 
field in Taurus at optical wavelengths (Moneti et al. 1984; 
Heyer et al. 1987; Goodman et al. 1990; Whittet et al. 
1992), but only a handful of polarization vectors exist 
in the near-infrared (Moneti et al. 1984; Tamura et al. 
1987; Goodman et al. 1992). Therefore, the magnetic 
field properties are poorly known in the high column den- 
sity regions where stars form. Near-infrared observations 
are needed to penetrate these regions. To address this 
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deficiency, in this paper we present 287 near-infrared po- 
larization measurements towards two regions of Taurus. 
One region is towards B213 and L1495 and has not been 
mapped before in polarization. Since L1495 is the highest 
column density area of Taurus with the highest density 
of embedded young stellar objects, measuring the mag- 
netic field in this region is important for studying the 
connection between magnetic fields and star formation. 
We also observed a diffuse region near where Goldsmith 
ct al. (2008) found striations in the 12 CO that matched 
the angles of the optical polarization vectors. 

In § 2 we discuss the observations, data reduction, and 
selection criteria to create the high-quality source cat- 
alogs. We also list references to previously published 
optical to near-infrared polarization studies in Taurus 
that will augment our new data set. We use these com- 
bined data sets to address several key magnetic field re- 
lated questions and issues in § 3. We discuss grain align- 
ment efficiency (§3.1), magnetic field strength (§3.2), 
cloud stability (§3.3), turbulence estimates in L1495 
(§3.4), and finally the magnetic field morphology in 
L1495 (§3.5). We summarize our findings in §4. 

2. OBSERVATIONS 

We used the Mimir instrument (Clemens et al. 2007) 
to observe -ff-band (1.6 /im) polarization of background 
starlight seen through the Taurus molecular cloud. 
Mimir is installed on the 1.8m Perkins telescope located 
near Flagstaff, AZ and operated by Lowell Observatory. 
The observations spanned the nights of 2009 Jan. 12- 
16 UT. We observed two regions within Taurus: one in 
a low-density portion of the cloud (hereafter 'Diffuse') 
and the other toward a high-density region, B213/L1495 
(hereafter 'Filament'). The Diffuse region is rectangular 
covering roughly 4 h 51 m — 4 h 53 m in Right Ascension and 
25°26' - 27°05' in Declination. The Filament region is 
not quite rectangular, but spans approximately 4 h 16 m to 
4 h 19 m in Right Ascension and 27°11' - 28°31' in Decli- 
nation. Coordinates are J2000. 

Each region consisted of an overlapping mosaic of 
10' x 10' Mimir fields-of-view. For each FOV, we ob- 
served in a six-position hex dither pattern where each 
dither position was offset by about 15" from the cen- 
ter. At each dither position we obtained images at 16 
angles of the half- wave plate (HWP), each separated by 
22.5°. The integration time per HWP angle was either 
2.5s (short) or 10s (long). The 16 HWP positions are 
equivalent to observing four angles separated by 22.5° 
four times each because the Stokes parameters vary as 
29 and because polarization "vectors" only have angle 
and not direction. Therefore, the effective integration 
time per pixel on the sky where all dithers overlap was 
4 x 6 x 2.5s = 60s (short) and 4 x 6 x 10s = 240s (long). 
These times are for each of the four independent HWP 
angles. The short integrations provide polarization mea- 
surements for most sources that were saturated in the 
long integrations. The Diffuse region has low extinction, 
therefore we observed it only with short integrations. For 
the Filament region, we observed all but five fields with 
both short and long integrations. The two regions ob- 
served are indicated in Figure 1 as hashed (short) and 
cross-hatched (short and long). All previously published 
optical to near- infrared wavelength polarization measure- 
ments are shown as white and black vectors. See § 2.2 for 



details on these catalogs. 

We reduced the data using custom IDL programs 
(Clemens et al. in prep) 5 . Each field, both long and short 
integrations, was processed separately to create stellar 
polarization lists. We then separately combined the short 
and long integrations from the different fields into a single 
catalog. For stars that position matched within 1" a star 
in adjacent field(s), we computed the weighted average 
of their Stokes parameters. Next we combined our short 
and long integrations in the Filament again using a 1" 
matching radius. For stars that matched we always used 
the polarizations obtained from the long integrations, as- 
suming they were more accurate. Stars saturated in the 
long integrations, however, did not have measured polar- 
izations. Therefore, when combining the short and long 
integrations, these stars have polarizations derived from 
the short integrations. 

2.1. High- Quality Vectors 

We applied two selection criteria to the catalogs to 
produce the high-quality polarizations used in this pa- 
per. First, we selected only those data with percentage 
polarization divided by its uncertainty p/a p > 3. Most of 
the sources removed by this criterion are faint and would 
also be removed by the second criterion. The low p/a p 
sources have median p w 3% and median a p « 2.5%. Sec- 
ondly, we required sources to be brighter than 12th mag. 
in the H band for the short integrations and brighter 
than 13th mag. (H band) for the long integrations. The 
second criterion eliminated the few fainter sources not 
caught by the first criterion. The measured polariza- 
tions for fainter sources are unreliable because they are 
dominated by errors not accounted for by the first cri- 
terion (D. Clemens, private comm.). These magnitudes 
are the 2MASS magnitudes for each source. Five sources 
in the Diffuse region did not have 2MASS counterparts 
and all five were excluded from the final catalogs. Nine 
of the stars in the Filament region did not have 2MASS 
counterparts and they were also excluded. 

The measured polarization from embedded protostars 
is not a useful probe of the magnetic field in molecular 
clouds. The radiation from the protostars may drive dust 
grain alignment, thus possibly altering the measured po- 
larization. Furthermore, reflected light from the disk and 
envelope is highly polarized and non-uniform. To avoid 
this potential bias, we compared our polarization posi- 
tions with those of the known protostellar members of 
Taurus from Luhman et al. (2006). None of our stars in 
the Diffuse region position matched within 1" a known 
Taurus member. These 125 sources are listed in Table 
1. In the Filament, 11 of the 173 polarization measure- 
ments matched within 1" the positions of known Taurus 
members. After excluding these, 162 high-quality po- 
larizations remained and these are listed in Table 2. Be- 
cause protostars may have some nebulosity leading to less 
certainty in their source positions, we also tried larger 
matching radii, up to 10", and found no difference in the 
number of protostars matched. 

Figures 2 and 3 show the vectors from the Diffuse and 
Filament regions overlaid on the 13 CO image from Gold- 
smith et al. (2008). The Mimir near-infrared polarization 

5 The latest released software can be found at: 
http: / / pcoplc.bu.edu/clemens/mimir / software.html 
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Figure 1. The two areas mapped are overlaid on the image of 3 CO emission (integrated over 2 — 9 km s~ 1 ) from Goldsmith ct al. 
(2008). Areas mapped with short integrations are shown with slanted lines; those mapped with both short and long integrations arc shown 
cross-hatched (§ 2). Also shown are polarization measurements from previous studies. 'Optical' are shown as white lines, 'infrared' as black 
lines, and 'I-band' are shown as thick black lines. Embedded sources in Taurus are shown as white circles (Luhman et al. 2006). Directions 
of increasing Galactic longitude (I) and latitude (b) are shown in the lower left. 



vectors are in black and the previously published optical 
vectors in white. In Figure 2 the Mimir vectors have po- 
larization angles that are similar to the optical ones. In 
Figure 3, some of the Mimir vectors are perpendicular to 
the B213 filament as are the optical vectors. Where the 
filament abruptly changes direction, the Mimir vectors 
also change direction and become oriented north-south. 
We include in Figure 3 a white dotted line separating the 
B213 filament (perpendicular vectors) and L1495 (paral- 
lel vectors). 

2.2. Previous Polarization Catalogs 

In §§3.2 and 3.3 we will combine previously published 
polarization measurements in Taurus along with the 
Mimir data. The published data arc identified by wave- 
length as follows: 'optical' (Moneti et al. 1984; Heyer et 
al. 1987; Goodman et al. 1990; Whittet et al. 1992), 'in- 
frared' (Moneti et al. 1984; Tamura et al. 1987; Goodman 
et al. 1992), and 'I-band' (Arce et al. 1998). The I-band 
designation is only approximate; the measurements were 
centered at 7660 A, with a bandpass of 2410 A. For the 
infrared, we only used if-band polarizations. Restrict- 
ing ourselves to K gave us the largest sample since some 



stars were only detected at K. We only used sources in 
these catalogs with p/cr p > 3. Finally, we excluded two 
vectors from Goodman et al. (1992) that are within 1" 
of known embedded sources (Luhman ct al. 2006), and 
two vectors from Arce et al. (1998) identified by them as 
being from stars closer than 150 pc, and thus likely in 
front of Taurus. 

3. RESULTS 

3.1. Alignment of Dust Grains in Dense Regions 

To examine how well the dust grains are aligned in 
different regions of Taurus, we plot in Figure 4 the per- 
cent polarization p versus column density, Ay, for the 
Mimir sources. To compute Ay, we obtained the 2MASS 
JHK S data for the Taurus cloud and used the NICER 
technique (Lombardi & Alves 2001) to create a 200" res- 
olution extinction map with 100" pixel spacing (Pineda 
ct al. 2010). NICER uses the 2MASS J-H and H-K s 
colors simultaneously to estimate the extinction towards 
each star. By using both colors, NICER has lower uncer- 
tainties than methods that use J — H or H — K s alone. 
The extinction map is then made by averaging together 
the Ay estimates for stars nearby each pixel. For each 
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Figure 2. Polarization map of the Diffuse region. Vectors from Table 1 are shown as black lines. See Fig. 1 for a description of the 
grayscale image and other vectors. For clarity, infrared vectors from previous studies arc not shown. 



polarization vector we use the Ay from the correspond- 
ing pixel in the extinction map. These values are listed 
in Tables 1 and 2. Estimating Ay from an extinction 
map produces more accurate values of Ay than just con- 
verting the measured color excess to Ay for each star 
because the unknown spectral type for each background 
star leads to large uncertainties in Ay. 

Following Goodman et al. (1992) and Whittet et al. 
(2008), we fit the data to a power law p = aA v b . Be- 
cause the errors in Ay are approximately all the same, 
we ignored them to find the best-fit least squares curves 
for the two regions, as shown in Figure 4. The dy- 
namic range of Ay in the Diffuse region is small enough 
that no correlation between polarization and Ay is seen 
(Pearson correlation coefficient r — 0.05). However, in 
the Filament region a significant relationship is found, 
p = (1.08 ± 0.06) V' 5 ^ - 04 , (r = 0.52). 

The three open circles in Fig. 4 arc not used in the 
fit because these have Ay values inconsistent with their 
measured color excess E(H — K). This can be seen in 
Figure 5 where we plot E(H — K) versus Ay for the two 
regions. To compute E(H — K) we assumed an average 
intrinsic color (H — K) = 0.114 ± 0.074 as measured by 
Pineda et al. (2010). The black lines show the expected 
extinction for a given E(H — K) assuming a standard 
dust extinction law (Rieke & Lebofsky 1985). In the Dif- 



fuse region, all the stars have E(H — K) within 2a of the 
black line. However, three sources in the Filament region 
are > 3c away from the black line. The low resolution 
of the extinction map has likely caused these sources to 
be assigned inaccurate Ay values. These three sources 
were thus excluded in the power law fit. 

Goodman et al. (1992, 1995) and Arce et al. (1998) 
measured the near-infrared and I-band polarization to- 
ward several regions in Taurus and Ophiuchus. They 
found that the percentage polarization had a slight de- 
pendence on Ay, but noted their data were also consis- 
tent with no dependence on Ay . They interpreted these 
results to argue that while increasing Ay added more 
dust grains that contribute to the column density, these 
grains did not add polarization. The lack of increased 
polarization with increasing Ay may be attributed to 
changes in dust properties (e.g., formation of icy mantles 
or grain growth), or lack of grain alignment in denser re- 
gions. They argued that measurements of the magnetic 
field by polarized absorption of background starlight only 
probe the magnetic field in a thin surface layer of the 
cloud. 

Whittet et al. (2008), however, used ~ 60 polarizations 
drawn from many datasets, including some with extinc- 
tions higher than those used by Goodman et al. (1992, 
1995) and Arce et al. (1998), and found a power law re- 
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Right Ascension (J2000) 

Figure 3. Polarization map of the Filament region. Vectors from Table 2 are shown as black lines. See Fig. 1 for a description of the 
grayscale image, vectors, and white circles. For clarity, infrared and I-band vectors from previous studies are not shown. The dashed white 
line denotes the division we made between the L1495 and B213 subregions (Fig. 6). 



lationship, p oc A V 0AS . They argued that the Radiative 
Torques mechanism (RT) for grain alignment (Dolginov 
& Mitrofanov 1976; Draine & Weingartner 1996, 1997; 
Lazarian & Hoang 2007) is the best explanation for their 
results. RT assumes there is an anisotropic radiation 
field impacting the dust grains. This radiation will im- 
part a net torque to the grains, assuming the grains are 
irregular in shape. For smooth grains with no helicity the 
torque will average out. The dust grains precess about 
the magnetic field lines and one component of the radia- 
tive torque is perpendicular to the rotation axis of the 
grain and points towards the center of the precession cir- 
cle, i.e. toward the magnetic field line. On average, this 
component of the torque acts to align the rotation axis of 
the grain with the magnetic field. Whittet et al. (2008) 
computed the efficiency for grain alignment by RT with 
a simple model. They assumed a spherical homogeneous 
cloud with an MRN dust grain size distribution (Mathis 
et al. 1977) and with an interstellar radiation field that 



matched the average from 0.1 — 100 /zm (Mathis et al. 
1983). They then computed the resultant polarization 
efficiency versus Ay and found that RT is capable of re- 
producing the observed data up to at least Ay = 10 mag. 
Beyond this critical Ay , p versus Ay should become flat 
unless grain growth occurs. 

Our power law index in the Filament region, 0.52 ± 
0.04, is consistent with the result obtained by Whittet et 
al.. Furthermore, our data are concentrated in a single 
small region while the data compiled by Whittet et al. 
(2008) span <~ 7° in Taurus. The agreement between our 
correlation and that found by Whittet et al. (2008) sug- 
gests that RT may be the primary mechanism for grain 
alignment in L1495. Furthermore, because we do not ob- 
serve a flattening of the p vs. Ay distribution, it implies 
the polarizing dust grains are not confined to a thin sur- 
face layer of the cloud, but that grain alignment occurs 
at least up to Ay of ~ 9 mag. This result is in contrast 
to Arce et al. (1998) who found a break point in p vs. Ay 
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Figure 4. Percent polarization versus Ay for: (a) the Diffuse region and (b) the Filament region. The solid curves are the least-squares 
fit to the data in each region, excluding the open circles in the Filament region. The fits are (a) p = (f .88 ± 0.09)Av — °'° °' 14 and (b) 
p = (f .08 ± 0.06)yly ' 52±0 - 04 .The black rectangle in (b) denotes the range plotted in (a). 
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Figure 5. E(H — K) versus Ay for: (a) the Diffuse region and (b) the Filament region. The black lines are the expected Ay for a given 
color excess, Ay = lt>.87E(H — K), assuming a standard extinction law (Rieke & Lebofsky 1985). The open circles in the Filament are 
the same sources as in Fig. 4. The black rectangle in (b) denotes the range plotted in (a). 



beyond which the distribution was flat at Ay = 1.3 ± 0.2 
mag. The difference in results is possibly explained be- 
cause the L1495 region has a higher column density than 
the region probed by Arce et al. (all their data have 
Ay < 4 mag). This is important because the efficiency 
of RT should increase with grain size and even more so 
when the wavelength of the radiation is comparable to 
the grain size (Cho & Lazarian 2005; Lazarian & Hoang 
2007). Grain growth is a collisional process, and so oc- 
curs more rapidly in regions with higher number density. 



Furthermore, in higher Ay regions, shorter wavelengths 
are extincted, but infrared photons may still align the 
large dust grains (Whittet et al. 2008). Lastly, RT re- 
lies on anisotropic radiation to work. Dense regions will 
be less fully penetrated by the interstellar radiation field 
than more diffuse regions, which may largely neutralize 
RT in the latter. 

Lastly, we should address the possible effect of proto- 
stars on our results. Whittet et al. found some evidence 
for increased polarization efficiency towards embedded 
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protostars. It was for this reason that we excluded vec- 
tors from stars that position-matched to known proto- 
stars in §2.1. However, the radiation from these pro- 
tostars could increase the measured polarization from 
nearby regions via RT. We consider a typical low-mass 
cloud core radius to be <~ 10,000 AU (e.g., Benson & 
Myers 1989). At the distance of Taurus 10,000 AU is 
~ 70". We then compared the polarization data in the 
Filament region to the embedded star catalog of Luhman 
ct al. (2006), using a 70" matching radius and found only 
two vectors match (with offset distances of 28" and 60"). 
Removing these two sources does not alter the best-fit 
relationship between p and Ay. Therefore, the increase 
in percent polarization with Ay is not an artifact caused 
by embedded protostellar illumination. 

3.2. The Magnetic Field Strength 

Although the total magnetic field strength, B, cannot 
be measured from our data, its plane-of-sky component, 
Bu , can be estimated from the dispersion in angle of the 
polarization vectors. Under conditions of flux-freezing, 
the magnetic field lines should be dragged inward as a 
region gravitationally contracts. Therefore, we expect B\\ 
should increase in higher density regions. For this reason, 
we divide the Taurus cloud into the subregions shown in 
Figure 6. These subregions conform with known regions 
in Taurus and the approximate boundaries of the 13 CO 
emission. The exact division between B213 and L1495 
can be seen more clearly in Figure 3, and corresponds 
to the Declination where the filamentary 13 CO emission 
changes direction abruptly. It is not exactly a straight 
line because we wanted to make it unambiguous which 
vectors were in each subregion. In addition to the cloud 
subregions, we chose two off-cloud subregions denoted 
OC1 and OC2. Furthermore, we also use the Diffuse and 
Filament Mimir data by themselves to make up a total of 
nine subregions. Note that the Diffuse subregion is just 
a subset of OC1. 

In this section we first present the derivation of the ve- 
locity dispersion and density of the gas in each subregion. 
These two quantities are needed to compute B\\. Then 
we will estimate Bu in each subregion using the Chan- 
drasekhar & Fermi (1953) and Hildcbrand ct al. (2009) 
methods. The results from this section are summarized 
in Tables 3 and 4. 

To estimate the velocity dispersion in each subregion, 
we first averaged the 12 CO and 13 CO spectrum from 
each 20" pixel in the subregion, using the data from 
Goldsmith et al. (2008). We then fit these averaged 
12 CO and 13 CO spectra with a gaussian to determine 
the radial velocity dispersion, <r(v). The line profiles in 
L1506 have two gaussian components separated by ap- 
proximately 3.5 km s _1 . Therefore, we fit two gaussians 
and average their dispersions. The individual gaussians 
have dispersions of 1.05 and 0.87 km s _1 for 12 CO and 
0.47 and 0.52 km s _1 for 13 CO. The uncertainties in ve- 
locity dispersions listed in Table 3 are from the fits. In 
this paper we use the 13 CO velocity dispersions to com- 
pute B\\, as 12 CO traces only a surface layer but 13 CO is 
sensitive to the whole cloud depth. As discussed in § 3.1, 
the polarizing dust does not appear confined to a thin 
surface layer. If we used the 12 CO velocity dispersion 
instead, all our B» estimates would increase by the ratio 



Table 3 

Average Physical Parameters for Subregions 



Velocity Dispersion 





A v 


n(H 2 ) 




linn 


Subregion 


(mag) 


(cm -3 ) 


(km s ) 


(km s ) 


B18 


1.80 ±0.10 


660 ± 50 


1.20 ±0.02 


0.90 ±0.02 


B213 


2.05 ±0.14 


750 ± 30 


1.16 ±0.02 


0.85 ±0.01 


Diffuse 21 


1.30 ±0.03 


200 ± 10 


0.84 ± 0.02 


0.61 ±0.02 


Filament a 


2.87 ±0.15 


920 ± 20 


1.30 ±0.03 


0.98 ±0.02 


HC12 


5.61 ±0.90 


970 ± 30 


1.14 ±0.02 


0.76 ±0.01 


L1495 


3.32 ±0.18 


990 ± 10 


1.19 ±0.02 


0.81 ±0.02 


L1506 b 


1.50 ±0.14 


620 ± 60 


0.87 ±0.01 


0.52 ±0.01 


OC1 


1.14 ±0.03 


210 ± 10 


0.93 ± 0.02 


0.71 ±0.01 


OC2 


0.91 ±0.09 


360 ± 50 


1.34 ±0.04 


0.83 ±0.03 


a Region only contains Mimir data 






b LI 506 has two velocity components 


separated by 


~3.5 km s- 1 



The velocity dispersion listed here is the average of the dispersion 
from each component. 

of the 12 CO/ 13 CO velocity dispersion. 

Obtaining the density is more difficult. Although the 
extinction map from Pineda et al. (2010) can be con- 
verted into column density, the number density remains 
uncertain since the line-of-sight cloud thicknesses are un- 
known. However, we can estimate the number densities 
from the CO data. Goldsmith et al. (2008) masked the 
pixels in the Taurus maps into three regions: maskO (nei- 
ther 12 CO nor 13 CO detected in individual pixels), maskl 
(only 12 CO detected), and mask2 (both 12 CO and 13 CO 
detected). Pineda et al. (2010) compared the gas and 
dust in Taurus, including effects such as depletion, tem- 
perature variations, and CO ices. From their results, we 
can use the following estimates for n(77 2 ): maskO (« 100 
cm~ 3 ), maskl (300 cm -3 ), and mask2 (> 10 3 cm" 3 ). 
The last is only a lower limit since n — 10 3 cm~ 3 is 
the critical density of 13 CO. For each polarization vec- 
tor, we then determine whether it lies in the maskO, 1, 
or 2 regions, and use the appropriate number density. 
Lastly, we compute the average n(H 2 ) for each subre- 
gion, using the standard deviation of the mean to repre- 
sent the uncertainty. HC12 and L1495, both sites with 
star-formation, show the highest average densities, close 
to 10 3 cm~ 3 , while the off-cloud subregions have lower 
densities of only a few hundred. 

An alternative, more geometric, approach to estimat- 
ing the density is to assume that the span on the sky of a 
subregion is approximately the same as the cloud depth 
in that subregion. We estimate maximum and minimum 
sizes for each subregion using the Ay = 2 mag. con- 
tour in the extinction map (Pineda et al. 2010). We use 
the Ay = 2 mag. contour because it outlines the high- 
density regions without picking up the low-level diffuse 
Ay. We then compute the geometric mean of the max- 
imum and minimum extents to arrive at an estimate of 
the cloud depth. If we assume a standard dust extinc- 
tion law (Rieke & Lebofsky 1985), N{H 2 ) = 9.4x 10 20 A V 
(Bohlin ct al. 1978). In all subregions, the geometric es- 
timates of n(H2) are less than those obtained from the 
masking method. For approximately square subregions 
like HC12 and L1495, the geometric densities are ~ 1/2 to 
1/3 of the mask-derived estimates, but for extended re- 
gions like B18 and B213, the geometric number densities 
are ~ 1/4 to 1/3 of the mask estimates. Furthermore, 
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Figure 6. Same as Fig. 1 with Mimir data shown as black vectors. For clarity, the embedded sources are not shown. We divided the map 
into the subregions labeled here, where HC12 is Heiles' Cloud 2, and OC1 and OC2 are two off-cloud regions. 



the geometric method cannot be applied to the two off- 
cloud regions and the Diffuse region. For these reasons, 
in this paper we use density estimates derived from the 
CO mask regions. 

3.2.1. The Chandrasekhar- Fermi Method 

The method of Chandrasekhar & Fermi (1953, here- 
after C-F) estimates the strength of the plane-of-thc- 
sky component of the magnetic field, B» , using the mea- 
sured dispersion in the angles of the polarization vectors. 
These authors assume the dispersion in angles arises from 
turbulence in the gas and that a strong magnetic field 
will resist the turbulence and hence will show a smaller 
dispersion in angle. They derived the expression: 



(1) 



where <j(&) is the unweighted dispersion in polarization 
angle, a(v) is the line-of-sight velocity dispersion of the 
gas, p is the gas density, and / is a factor that accounts 
for averaging of the magnetic field direction along the 
line-of-sight. The last is important because along any 
given line-of-sight there will be many turbulent cells each 
with a different plane-of-sky polarization angle. The 
measured polarization angle at any point is then the 



average polarization along the line-of-sight through the 
cloud. This line-of-sight averaging will lead to a smaller 
measured a (9), and thus result in an overprediction of 
B«. 

Using Equation 1 and the values in Table 3, we es- 
timate Bii in each subregion. To compute p we con- 
vert the number densities in Table 3 via the equation 
p = n(H2)mHp,H 2 , where tuh is the mass of the hydro- 
gen atom and pn 2 ~ 2.8 is the mean molecular weight 
per hydrogen molecule. Lastly, we set / = 0.5. Ostriker 
et al. (2001) found from simulations that this value for 
/ is a reasonable approximation provided <r(0) < 25°. 

Uncertainties in B^ are computed from propagation 
of errors on Equation 1. The propagation of errors is 
straightforward except for the uncertainty in a (6), which 
is itself a measured dispersion. We go back to the def- 
inition of standard deviation and propagate the errors 
to obtain the following expression for the uncertainty in 

a(ey. 



1 



(N-l)a(0) 



(2) 



where 9i and a$ i are the angle and uncertainty for each 
polarization vector and is the average angle for all the 
vectors in a subregion. 
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Table 4 lists the values for Bu in each subre- 
gion. The magnetic field strengths are 10 — 17 fj,G in 
the low-density off-cloud subregions, ~ 25 /iG in the 
B213/Filament/L1495 subregions, and peak at 42±4/iG 
in HC12. These values are consistent with previous 
estimates for magnetic field strengths. Heyer et al. 
(2008) used data from part of OC1 and estimated Bn = 
14 /iG based on magnetohydrodynamic induced velocity 
anisotropy and separately also from the C-F method. 
This value is in excellent agreement with our values of 
17 /iG and 12 /iG for the diffuse and OC1 regions, re- 
spectively. The line-of-sight component of the magnetic 
field, Bi os , has also been measured using OH Zeeman 
line splitting. Troland & Crutcher (2008) measured OH 
Zeeman line splitting towards 34 cores, including 11 in 
Taurus, but only 2 of these 11 with a significance > 3a. 
Those two are B217-2 (B los = 13.5i3.7juG) and TMC1 
{Bi os = 9.1 ± 2.2 /xG). TMC1 is a core within our HC12 
subregion. The much smaller Bi os measurement com- 
pared to Bn implies that the magnetic field lies close to 
the plane-of-the-sky. However, the relatively small num- 
ber of polarization vectors in HC12 (22) means there may 
be systematic errors in the measured value of B\\ . 

L1506 and OC2 have a(9) much larger than 25°, mean- 
ing that our assumption of / = 0.5 may be inaccurate, 
and thus B^ is less certain for these subregions. L1506 
stands out with a(9) = 63° and the smallest B\\. In Fig- 
ure 6 the vectors in the western half of L1506 are approx- 
imately parallel to the cloud while those in the eastern 
half are more perpendicular, yielding a large a(9) and 
correspondingly small B\\. Future studies should treat 
these two regions separately. 

3.2.2. The Hildebrand et al. Method 

The C-F method has the advantage of being straight- 
forward to implement, but it will yield lower limits for 
.B|| because it assumes that the dispersion in polarization 
angle is entirely due to turbulence. Large-scale, non- 
turbulent changes in the magnetic field direction need to 
be accounted for because Taurus has large-scale compo- 
nents (see e.g., Fig. 1). The method of Hildcbrand et 
al. (2009) accounts for non-turbulent variations without 
assuming a model field. This method starts by consid- 
ering a two-dimensional map of the magnetic field pro- 
jected on the plane-of-the-sky, where at any position x, 
the angle of the magnetic field is $(x). The difference 
in angle $(x) — $(x + £) is then computed for every pair 
of vectors. These differences in angle are then binned 
by distance, £, and the sum over the N{£) pairs of vec- 
tors for that bin is computed to arrive at the two-point 
correlation (called the dispersion function): 



(A$ 2 (£)) 1/2 = 



N(i) 



£[$(x)-$(x + *)]2. (3) 



Hildebrand et al. assume that B(x) is composed of 
a large-scale structured field, B (x), and a turbulent 
component, Bt(x). Because Bo(x) is a smoothly vary- 
ing quantity, its contribution to the dispersion function 
should increase linearly with I for small distances. The 
turbulent component, Bb(x), should also increase with 
I up to a maximum value when £ exceeds the turbulence 



correlation length S. Through Taylor series expansion, 
Hildebrand et al. separate the contributions from Bo (x) 
and B t (x) to the dispersion function. The authors show 
that the square of the dispersion function can be approx- 
imated as: 



(A<S> 2 (£)) tot = b 2 +m 2 £ 2 +a 2 M (£), 



(4) 



where (A$ 2 (f )) to t is the dispersion function computed 
from the data. The quantity <j 2 M {£) is computed from the 
A$ (£) = $ (x) - $ (x + £) values in each bin. Each A$(£) 
has an associated variance obtained from propagation of 
errors. Because the dispersion function is computed by 
summing the square of the difference in angles, measure- 
ment errors are added in quadrature and thus will bias 
the computed (A<fr 2 (£)} to t by u 2 M {£), where o 2 M {£) is sim- 
ply the average of the variances on A$(£) in a bin. 

The quantity b 2 is the intercept of a straight line fit to 
the data (after subtracting As the distance, £, 

approaches zero the contribution to (A$ 2 (£)) from B (x) 
disappears (represented by m 2 ) and only B t (x) remains 
(represented by b 2 ). Equation 4 is valid for displacements 
£ larger than the correlation length for B t (x), S, and for 
£ much smaller than the length scale for variations in 
B (x), d. 

Hildebrand et al. then solved their equation for b 2 to 
find the ratio of the turbulent to the large-scale magnetic 
field strength: 



(jB?) i/ 2 



B n 



VT^b 2 



(5) 



Lastly, Hildebrand et al. assumed that a{9) ~ SB/B , 
where SB is the variation in magnetic field about the 
large-scale field Bg. Therefore, making the inference that 
{B 2 ) 1 / 2 corresponds to SB, Equation 5 is an expression 
for a(theta). 

We implement the method of Hildebrand et al. (2009) 
as follows. We first compute the difference in polariza- 
tion angle between every set of two points in a subregion. 
Next, we bin these data into either 3' or 5' wide bins and 
compute the average dispersion function in each bin. We 
chose bin widths of 3' or 5' to have N(£) > 10 for the 
bins that are fit with a straight line. Only one bin has 
< 10 measurements, this is the third bin in HC12 with 
N{£) = 8. The five subregions with Mimir data (B213, 
Diffuse, Filament, L1495, and OC1) have many more 
points than the other four, allowing us to use a 3' bin 
size for these subregions. Next, we subtract the average 
variance, o 2 M {£), in each bin. Lastly, we fit a straight 
line to the data versus distance squared. We always ex- 
clude the first bin because in all subregions this bin has 
~ 1/3 the number of points compared with the remain- 
ing bins. Furthermore, excluding the first bin ensures 
that our distances will be greater than S. Houde et al. 
(2009) estimated S = 0.016 pc towards OMC-1. At the 
distance of Taurus 0.016 pc is 24". With the exception 
of three subregions, discussed in the next paragraph, the 
maximum distance fit was either 20' or 21', depending 
on the bin width. This maximum distance minimizes 
the chance that £ > d. The zero distance intercept of the 
straight line fit is b 2 in Equation 4. Figure 7 shows our 
results from using the Hildebrand et al. method. The 
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Table 4 

Magnetic Parameters for Subrcgions 



B]± £ 

b 2 # of a(6) C-F Hildcbrand C-F Hildcbrand 

Subrcgion (deg 2 ) vectors (deg) (pG) 



B18 


527 ± 52 


62 


26 


19 ± 1 


30 ± 2 


0.22 ±0.02 


0.14 ±0.01 


B213 


192 ± 24 


142 


18 


28 ± 1 


51 ± 4 


0.17 ±0.01 


0.09 ±0.01 


Diffuse a 


75 ± 18 


125 


11 


17 ± 1 


31 ± 4 


0.18 ±0.01 


0.10 ±0.01 


Filament a 


124 ± 32 


162 


24 


27 ± 1 


82 ± 11 


0.25 ±0.02 


0.08 ±0.01 


HC12 


144 ± 80 


22 


13 


42 ± 4 


61 ± 17 


0.32 ±0.06 


0.22 ±0.07 


L1495 


103 ± 32 


99 


23 


25 ± 1 


77 ± 12 


0.32 ±0.02 


0.10 ±0.02 


LI 506 


921 ± 4 


39 


63 


5 ± 1 


12 ± 1 


0.79 ± 0.08 


0.29 ±0.03 


OC1 


73 ± 23 


191 


18 


12 ± 1 


37 ± 6 


0.22 ±0.01 


0.07 ± 0.01 


OC2 


158 ± 38 


40 


34 


10 ± 1 


38 ± 6 


0.21 ±0.03 


0.06 ±0.01 



a Region only contains Mimir data 



bins are plotted versus distance (not distance squared) 
and the uncertainty in each bin is the unweighted stan- 
dard deviation of the mean. The black line shows the 
best fit and its length represents the range of distances 
used in the fit. 

Our expectation is that (A<J> 2 (£)) should increase ver- 
sus offset distance as Bt(x) and Bo(x) will cause a larger 
dispersion in angles at larger separation distances (see 
also Fig. 1 in Hildcbrand et al. 2009). However, in three 
subregions the best-fit line had a negative slope when 
fitting over the range 5' - 20' (B18) or 3' - 21' (Diffuse 
and OC1). For these three subregions we increased the 
range of distances used in fitting until the slope was pos- 
itive at the > lcr level. For B18 the final distances fit 
are 5' — 35' and for the Diffuse and OC1 subregions the 
range is 3' — 36'. These are the ranges shown in Figure 7, 
and the b 2 values in Table 4 are derived from them. By 
using the larger distances in fitting we run a greater risk 
of I exceeding d. The approach we use may be justified 
in the Diffuse and OC1 subregions, where from Figure 1 
it appears that d is greater than 1°, but it is unclear if 
the same is true for B18. However, the important quan- 
tity for our purposes is the intercept, 6 2 , which is less 
sensitive to the maximum fit distance used, giving us 
confidence in our best-fit values. In B18, the best-fit b 2 
varies by less than the error for maximum fit distances 
25' — 45', and in the Diffuse and OC1 subregions, maxi- 
mum distances from 30' to 60' produce best-fit values of 
b 2 that vary from each other by less than the error listed 
in Table 4. 

With our b 2 values, we use Equation 5 to obtain <j{9) 
and finally Equation 1 to estimate Bu . Our results are 
listed in Table 4. The Hildebrand et al. method pro- 
duces estimates of Bu that are <~ 1.5 — 4x larger than 
those from C-F. This increase is expected as the former 
removes large-scale B-ficld effects. In general, denser re- 
gions tend to have higher values for Bu . This behavior is 
not surprising since Bu oc y/p from Equation 1, but may 
also indicate that frozen magnetic field lines are being 
dragged inward by gravity. 

The relatively limited number of vectors in regions 
without Mimir data may introduce biases in our esti- 
mates of For example, L1506 has an extremely high 

value for b 2 that is determined to high precision, lead- 
ing to a well-determined estimate for Bu that is much 
lower than what is found even in low-density regions 
(12 ± 1 [iG). Another possible source of systematic er- 



ror is the correction factor /. Because Hildebrand et 
al. (2009) accounts for non-turbulent variations in mag- 
netic field direction, it is possible that / w 0.5 even in 
regions where a{9) > 25°. It is also possible that a dif- 
ferent correction factor is needed for the Hildcbrand et 
al. technique compared to C-F. More work is needed to 
understand how / depends on physical quantities. 

3.3. Cloud Stability 

Is the magnetic field strength in Taurus sufficiently 
strong to be important in regulating star formation? 
Nakano & Nakamura (1978) derived the stability crite- 
rion for an isothermal gaseous layer threaded by a per- 
pendicular magnetic field to be (M/$) cr it = l/V47r 2 G. 
For a mass-to-magnctic flux ratio exceeding this value, 
the region will collapse owing to gravity. The stability 
criterion can be rewritten in terms of the dimensionless 
magnetic critical index: 

(M/$) 

^wm^r 7 - 6N ^ )/B ^ (6) 

where Nn (#2) is the column density in units of 10 21 cm~ 2 
along a magnetic flux tube and B tot is the total magnetic 
field strength in /iG. For fi > 1, the flux tube is super- 
critical, meaning it should collapse due to self-gravity, 
but for /i < 1 the tube is magnetically supported. 

Because of projection effects between N\\ (H2)/ B to t and 
the observed quantity N{H2)/B\\, fj, b s will overestimate 
fj, with an average correction factor \x = /U 6s/3, assuming 
a random orientation of the magnetic field with respect to 
the line of sight (Heiles & Crutcher 2005). Furthermore, 
for our data we have measurements of Ay, not N{Hi). 
The conversion factor we used before is N(H 2 ) = 9.4 x 
10 20 Ay. Putting these two together, Equation 6 can be 
rewritten as: 

fi = 2AA v /B h (7) 

with Ay in magnitudes and B11 in /jG. 

Using the Ay values from Tables 1 and 2, we compute 
the weighted average Ay and the weighted standard de- 
viation of the mean in each subregion. The standard 
deviation of the mean is larger than the statistical uncer- 
tainty in all subregions. These average values are listed 
in Table 3. Using Equation 7, we compute \x for each 
subregion with both the C-F and Hildebrand et al. esti- 
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Figure 7. Plot of the dispersion function versus distance for the nine subregions. The black line shows the best-fit power law to the data 
points. Its length denotes the range of distances used in fitting. We subtracted <J%i(£) from (A«I> 2 ) before plotting. 
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However, the connection between the two is important. 
The magnetic critical index only regulates collapse due 
to ambipolar diffusion. It is still possible for gas to col- 
lapse along the field lines. The B213 filament appears to 
be a prime candidate for this scenario since the magnetic 
field is parallel to the short axis of the filament along 
its entire length. Small cores within the filament, not 
probed by our observations, would then be able to form 
stars. 



Figure 8. Plot of the magnetic critical index, fi, from Table 4 for 
each subregion. The weighted average values (excluding L1506) for 
the C-F and Hildebrand et al. techniques are plotted as black and 
gray dashed lines, respectively. 

mates of B\\ . The results are listed in Table 4 and plotted 
in Figure 8. 

All subregions are subcritical. The average /i, exclud- 
ing L1506, is 0.21 (C-F) or 0.09 (Hildebrand et al). It is 
remarkable that \i is roughly constant across subregions. 
This constancy implies the bulk of the cloud is not in- 
volved with star formation. In § 3.1 we argued that our 
polarization measurements probed up to at least a col- 
umn density Ay <~ 9 mag. Since even subregions with 
known protostars are highly subcritical, it appears our 
measured polarizations and estimates of B\\ may not be 
probing the innermost regions where stars are forming. 



3.4. Turbulence 

In § 3.1 we argued that our observed power law in the 
Filament region, p oc Ay ' 52 ^' 04, , was Qi ue to the Radia- 
tive Torques model for grain alignment. Whittet et al. 
(2008) did some simple modeling to show that RT was 
effective even at high Ay. However, they did not con- 
sider the possible effects of turbulence, which may also 
explain this behavior. Jones et al. (1992, hereafter JKD) 
combined data from various lines-of-sight to obtain po- 
larization at K (2.2 ^m) versus optical depth tk- Using 
the relation tk = 0.09Ay from JKD, we fit their data 
to a power law to obtain px oc A v °- 79 , which is steeper 
than our observed power law. JKD did not model the 
grain alignment mechanism, but did include the effects 
of turbulence with two models. Their wave model consid- 
ered a superposition of cloudlets along the line-of-sight, 
each threaded by an Alven wave with random phase and 
plane of vibration. All waves had the same direction of 
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propagation. Their component model was an extension 
of the model of Myers & Goodman (1991). This model 
assumed a magnetic held with a constant uniform compo- 
nent plus a random component that could have arbitrary 
angles for different turbulent cells along the line-of-sight. 
JKD found that both models best fit their data when 
cquipartition existed between the magnetic and turbu- 
lent energy densities (wave model) or when the strengths 
of the uniform and random components of the magnetic 
field are equal (component model). Both of their mod- 
els predict that we should measure stronger turbulence 
because our power law index is lower than theirs. 

We test the predictions of JKD by computing the rel- 
evant ratios for the wave and component models in the 
Filament region from the values listed in Tables 3 and 
4. For the wave model, the turbulent to magnetic energy 
density ratio can be written as the ratio of the turbulent 
to Alven velocity, ct(v)/Va, where Va = B/y/iirp. Using 
the 1 CO velocity dispersion and Bu from the Hildebrand 
et al. method, we find that a(v)/V A = 0.39 ± 0.04. For 
the component model, the ratio of the turbulent-to-large- 
scale field strength, (B 2 ) 1 / 2 / B , is exactly computed by 
the Hildebrand ct al. method via Equation 5. In the Fil- 
ament, (B 2 y/ 2 /B = 0.20 ±0.02. Therefore, both ratios 
are significantly less than 1 and turbulence is less im- 
portant than magnetic fields. This finding contradicts 
the predictions of the JKD models. Both JKD models 
assume a constant magnetic field component that does 
not vary in angle, and all dispersion in angle is due to 
turbulence. However, we know this is not the case in Tau- 
rus. The Filament region analyzed has a broad curved 
shape that does not appear to be caused by turbulence. 
In the next section we examine possible causes for this 
morphology. 

3.5. Magnetic Field Morphology 

Figure 3 shows the B213 filament and L1495 region 
with polarization vectors overlaid. The filament is re- 
vealed by the 13 CO emission, which serves as a dense 
gas tracer. In B213, the magnetic field is perpendicular 
to the apparent long axis of the filament. This mor- 
phology can be readily explained if the gas and dust 
have gravitationally collapsed along the field lines to 
form the filament. Between B213 and L1495, the fila- 
ment turns sharply and becomes oriented approximately 
north-south. The magnetic field also abruptly transitions 
from being perpendicular to being parallel to the fila- 
ment. Furthermore, the filament itself appears slightly 
curved, a curvature that is well-matched by the magnetic 
field. 

The behavior of the magnetic field in L1495 can be ex- 
plained if this region is the compression front of a bubble. 
It would explain the curved nature and the fact that the 
magnetic field and gas are parallel to each other. Un- 
der conditions of flux-freezing, an expanding bubble will 
carry the magnetic field with it and stretch it to be par- 
allel to the edge of the compression front (e.g., Novak et 
al. 2000; Li et al. 2006). 

In Figure 9 we show 12 CO, 13 CO, and H I emission in 
L1495 integrated from 7.6 — 8.9 km s _1 . This velocity 
range shows the curved bubble-like nature of this region 
more clearly than the 2 — 9 km s _1 range used in the 
previous figures. The 12 CO and 13 CO are from Gold- 



smith et al. (2008) and the H I is Arecibo data from M. 
Krco (PhD Thesis, Cornell University, in preparation). 
In each panel, a curved structure in the gas is present 
that matches the location and radius of curvature of the 
polarization vectors. Note that the B213 filament is not 
evident in the 12 CO and 13 CO emission because it emits 
at lower local standard of rest (LSR) velocities. 

Figure 10 shows the shell of a bubble-like structure 
visible in the H I emission to the northwest of L1495, 
outside the area mapped in 12 CO and 13 CO. The H I is 
integrated from 8.5 — 9.6 km s" 1 to show both the near 
and far edges of the shell. This shell is approximately 7° 
in diameter, with one edge at L1495 and the other edge 
at the top of the map. The H I gas at the top edge of the 
shell is more redshifted than the gas in L1495, implying 
that whatever source has created the shell is located at 
a distance comparable to or greater than that of Taurus. 

O and B stars may have sufficient stellar fluxes or 
winds to clear the interstellar medium and create the 
observed bubble. Using SIMBAD 6 we select all O and B 
stars within the area of Figure 10 that also have paral- 
lax angles of 3 — 7 milliarcseconds (parallax distances of 
142 - 333 pc). No O stars are in the field, and only 8 B 
stars. These are shown as squares in Figure 10. None of 
the B stars are likely candidates for creating the bubble. 
Their positions are scattered throughout the map, and 
they are too distant. Only one B star is closer than 200 
pc (at 180 pc) compared with the ~ 140 pc distance of 
Taurus. This star is at 3 h 46 m +29°40', but the closest B 
star in angular separation to L1495 (4 h 03 m +28°07') is 
at 300 pc. 

Alternatively, the shell could be a supernova rem- 
nant. A recently discovered gamma-ray pulsar, Fermi- 
LAT PSR J0357+32 has been found in the field by Abdo 
et al. (2009). Using the improved source position from 
Ray et al. (2011), we marked the position of this source 
with a cross in Figure 10. This source has been sought in 
radio with both Arecibo and GBT, but not detected (Ray 
et al. 2011). Current models of pulsars have a gamma- 
ray beam that is broader than the radio beam (Watters 
et al. 2009), so Earth may simply not be within the ra- 
dio beam cone. Because this source has only recently 
been discovered, little is known about it. The improved 
positional accuracy of Ray et al. (2011) leads to more 
accurate estimates of the spin-down frequency ,z>, and a 
characteristic age r = —v/i> — 540 kyr. As a pulsar that 
is possibly nearby, with a candidate supernova remnant 
interacting with the Taurus molecular cloud, it should 
be studied further. 

4. SUMMARY AND CONCLUSIONS 

We have observed near-infrared .ff-band polarization 
of background starlight seen through two regions of the 
Taurus molecular cloud using the Mimir instrument, 
mounted on the 1.8m Perkins telescope located near 
Flagstaff, AZ and operated by Lowell Observatory. After 
data reduction and selection of high-quality vectors, we 
obtained 125 vectors in the diffuse region and 162 vec- 
tors in the B213/L1495 region. The latter region had not 
been previously studied in polarization. Our results are 
as follows. 

6 http://simbad.u-strasbg.fr/simbad/ 
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Figure 9. Polarization map in the Filament region overlaid on different background images. (Left) and (middle) - 12 CO and 13 CO 
emission, respectively, from Goldsmith ct al. (2008). (right) -Hi from M. Krco (private comm.). All three emission images were integrated 
from ~ 7.6 — 8.9 km s _1 . The vectors and white circles are the same as in Fig. 3. 




04 h 24 m 18 m 12 m gin 

Right Ascension (J2000) 

Figure 10. H I emission integrated from 8.5 — 9.6 km s _1 with vectors overlaid as in Fig. 3. The squares denote the positions of B 
stars in the field having parallax distances between 142 — 333 parsccs, while the cross is the position of gamma-ray pulsar Fermi-LAT PSR 
J0357+32. 



Grain Alignment - In the B213/L1495 filament we 
found an increase in the percentage polarization, p, ver- 
sus column density, Ay. These data are well-fit by a 
power law with p oc A V °- 52±0M . The Radiative Torques 
model for dust grain alignment is able to explain the ob- 
served trend. Our near-infrared data probe the magnetic 
field geometry up to at least Ay ~ 9 mag., not just in a 
thin surface layer. 

Magnetic Field Strength - We divided our data, along 
with the ~ 400 previously published optical to near- 



infrared polarization vectors, into nine subregions within 
Taurus. In each subregion we estimated the strength 
of the plane-of-the-sky component of the magnetic held, 
£>l|, from the dispersion in polarization angles. The 
Chandrasekhar- Fermi technique (C-F) uses these disper- 
sions directly to estimate B\\ while the Hildcbrand et al. 
method accounts for large-scale spatial variations in po- 
larization angle not arising from turbulence. In all sub- 
regions we find the Hildcbrand et al. method produces 
larger estimates for B\\ than the C-F method. The mag- 
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nctic field in L1506 appears to have very different angles 
in the eastern and western halves leading to very small 
values for B^ . The Hildcbrand et al. method accounts for 
this variation better than C-F. Even so, the small num- 
ber of vectors in L1506 leads to an estimate for B\\ that 
appears anomalously small (12 ^.G) compared to other 
subregions. 

Cloud Stability - We used the derived estimates for B\\ 
to compute the critical mass to magnetic flux ratio index, 
fj,. In all subregions this index is < 1, indicating that the 
clouds are supported by the magnetic field. Furthermore, 
(i is nearly uniform across all the subregions with an 
average value of 0.21 (C-F method) or 0.09 (Hildcbrand 
et al. method). The constancy of /i even in high column 
density regions like L1495 further supports the idea that 
our near-infrared data may not be probing the magnetic 
field in the densest star-forming cores. However, \i only 
governs gravitational collapse across the field lines. The 
B213 filament is a prime example where it appears the 
gas and dust have collapsed along the field lines to form 
a filament. 

Cloud Turbulence - The power law fit to our data in 
B213/L1495 is flatter than the observed trend seen by 
Jones et al. (1992). Based on their models for turbu- 
lence, we expected to find a turbulent to magnetic field 
energy density ratio, cf(v)/Va, greater than 1, and also 
the ratio of the turbulent to large-scale magnetic field 
strength > 1. We calculated both ratios to be < 1 for 
Taurus, implying that turbulence is less important than 
the magnetic field. The component model of Jones et al. 
(1992) does not allow the non-turbulent component of 
the magnetic field to vary in angle. Furthermore, their 
wave model does not allow the Alven wave to vary in 
direction of propagation. The Hildebrand et al. method 
allows for non-turbulent variations in angle, and such 
behavior is observed in Taurus (see e.g. Figure 1). Non- 
turbulent variations in angle of the magnetic field may 
explain the difference between the predictions of Jones 
et al. (1992) and the calculated quantities. 

Magnetic Field Morphology - In the L1495 region the 
Mimir data show that the magnetic field appears to have 
a large-scale curvature. If the Mimir data arc tracing 
the shell of a bubble, its curvature would explain the 
observed morphology. Arc-shaped segments of a shell 
are seen in H I. We looked for potential causes of this 
structure but found no O stars exist in the region at the 
appropriate distance and only 8 B stars, none of which 
seems a likely candidate based on location and parallax 
distance. However, the shell may be a supernova rem- 
nant. A recently discovered gamma-ray pulsar, Fcrmi- 
LAT PSR J0357+32, is in the field. It has been searched 
for, but remains undetected at radio frequencies. How- 
ever, this non-detection may be a consequence of the 
radio beam being narrower than the gamma-ray beam. 

The Taurus molecular cloud remains an excellent tar- 
get for understanding the importance of magnetic fields 
in molecular clouds. Near-infrared polarimeters such as 
Mimir can map the large-scale magnetic field to moder- 
ate optical depth. In conjunction with upcoming submil- 
limeter polarimeters like ALMA and SCUBAPOL2, we 
will be able to probe the magnetic field direction from 
the largest to smallest spatial scales. 
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± 


2 


1.29 


± 


0.22 


88 


04 


51 


57.7 


26 


40 


28 


2.85 


± 


0.46 


60 


± 


5 


1.26 


± 


0.23 


89 


04 


51 


58.2 


27 


02 


29 


1.56 


± 


0.31 


59 


± 


6 


1.65 


± 


0.20 


90 


04 


51 


59.4 


25 


46 


08 


2.19 


± 


0.35 


22 


± 


5 


1.65 


± 


0.25 


91 


04 


51 


59.9 


26 


41 


22 


2.40 


± 


0.35 


52 


± 


4 


1.26 


± 


0.23 


92 


04 


51 


59.9 


25 


59 


06 


2.20 


± 


0.48 


38 


± 


6 


1.48 


± 


0.24 


93 


04 


52 


00.5 


26 


38 


36 


2.04 


± 


0.41 


70 


± 


6 


1.29 


± 


0.22 


94 


04 


52 


01.3 


26 


48 


59 


1.81 


± 


0.39 


44 


± 


6 


0.90 


± 


0.23 


95 


04 


52 


02.3 


26 


57 


28 


2.01 


± 


0.19 


53 


± 


3 


1.34 


± 


0.24 


96 


04 


52 


02.3 


26 


54 


53 


2.12 


± 


0.68 


45 


± 


9 


1.04 


± 


0.21 


97 


04 


52 


04.2 


25 


44 


11 


1.29 


± 


0.20 


22 


± 


4 


1.73 


± 


0.22 


98 


04 


52 


05.4 


26 


01 


54 


1.72 


± 


0.27 


47 


± 


4 


1.56 


± 


0.20 


99 


04 


52 


06.1 


25 


32 


00 


0.84 


± 


0.25 


58 


± 


9 


1.14 


± 


0.22 


100 


04 


52 


06.7 


25 


37 


12 


0.75 


± 


0.21 


58 


± 


8 


0.97 


± 


0.25 


101 


04 


52 


06.8 


26 


24 


56 


1.39 


± 


0.43 


48 


± 


9 


0.78 


± 


0.23 


102 


04 


52 


07.1 


26 


18 


42 


1.28 


± 


0.17 


54 


± 


4 


1.32 


± 


0.22 


103 


04 


52 


07.1 


25 


33 


48 


0.65 


± 


0.14 


47 


± 


6 


1.16 


± 


0.22 


104 


04 


52 


07.7 


26 


02 


30 


1.94 


± 


0.22 


38 


± 


3 


1.56 


± 


0.20 


105 


04 


52 


08.3 


26 


13 


52 


1.49 


± 


0.23 


46 


± 


4 


1.34 


± 


0.25 


106 


04 


52 


08.7 


25 


43 


36 


1.20 


± 


0.07 


34 


± 


2 


1.63 


± 


0.24 


107 


04 


52 


09.5 


25 


39 


50 


0.78 


± 


0.09 


40 


± 


3 


1.17 


± 


0.25 


108 


04 


52 


09.8 


26 


23 


03 


1.36 


± 


0.12 


55 


± 


3 


1.10 


± 


0.22 


109 


04 


52 


11.2 


27 


02 


15 


2.22 


± 


0.32 


56 


± 


4 


1.56 


± 


0.22 


110 


04 


52 


12.1 


26 


03 


39 


1.34 


± 


0.11 


41 


± 


2 


1.59 


± 


0.21 


111 


04 


52 


12.3 


26 


52 


06 


2.20 


± 


0.67 


61 


± 


9 


0.90 


± 


0.22 


112 


04 


52 


12.8 


26 


07 


35 


1.03 


± 


0.25 


40 


± 


7 


1.48 


± 


0.24 


113 


04 


52 


12.9 


26 


35 


53 


1.92 


± 


0.20 


58 


± 


3 


1.05 


± 


0.23 


114 


04 


52 


13.1 


26 


25 


02 


1.76 


± 


0.51 


53 


± 


8 


0.72 


± 


0.22 


115 


04 


52 


13.6 


26 


31 


21 


1.76 


± 


0.17 


61 


± 


3 


0.92 


± 


0.23 


116 


04 


52 


13.9 


27 


00 


38 


2.28 


± 


0.44 


52 


± 


6 


1.41 


± 


0.26 


117 


04 


52 


14.6 


25 


55 


24 


0.99 


± 


0.27 


43 


± 


8 


2.21 


± 


0.25 


118 


04 


52 


15.8 


26 


20 


53 


1.13 


± 


0.15 


68 


± 


4 


1.23 


± 


0.20 


119 


04 


52 


16.3 


26 


22 


03 


1.40 


± 


0.24 


66 


± 


5 


1.05 


± 


0.20 


120 


04 


52 


16.6 


26 


54 


04 


2.26 


± 


0.25 


50 


± 


3 


0.97 


± 


0.23 


121 


04 


52 


17.1 


26 


39 


18 


1.39 


± 


0.40 


62 


± 


8 


0.82 


± 


0.26 


122 


04 


52 


17.1 


26 


44 


01 


1.55 


± 


0.47 


58 


± 


9 


0.79 


± 


0.24 


123 


04 


52 


17.8 


26 


51 


26 


2.02 


± 


0.44 


53 


± 


6 


0.90 


± 


0.22 


124 


04 


52 


20.7 


25 


53 


03 


2.81 


± 


0.78 


29 


± 


8 


2.57 


± 


0.28 


125 


04 


52 


22.8 


26 


55 


59 


1.47 


± 


0.41 


62 


± 


8 


0.96 


± 


0.23 



Note. — Right Ascension, a, is given as hours, minutes, and seconds and Declination, 5, is given as degrees, minutes, seconds. 
a Angles are equatorial, measured east from north, 
k Extinction map from Pineda et al. (2010) 
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Table 2 

//-band Polarization Data for Stars in the Filament Region 



btar 




a 




S 








P 




6 


la 




A 


b 

V 


# 


(J2000) 


(J2000) 




(%) 


(D 


eg- 


) 


(mag 


;) 


1 


04 


16 


36.5 


27 42 


59 


3 


06 


± 


0.89 


79 


± 


8 


1.17 


± 


0.37 


2 


04 


16 


37.9 


28 16 


06 





.92 


± 


0.27 


89 


± 


8 


4.62 


± 


0.35 


3 


04 


16 


38.1 


27 52 


52 


1 


34 


± 


0.25 


45 


± 


5 


2.26 


± 


0.36 


4 


04 


16 


38.3 


27 36 


09 


1 


69 


± 


0.47 


53 


± 


8 


1.38 


± 


0.25 


5 


04 


16 


40.2 


27 37 


26 


6 


12 


± 


1.26 


86 


± 


6 


0.90 


± 


0.27 


6 


04 


16 


40.4 


28 29 


18 


1 


.19 


± 


0.22 


-11 


± 


5 


2.59 


± 


0.33 


7 


04 


16 


40.4 


28 26 


59 


1 


03 


± 


0.20 


-7 


± 


6 


3.06 


± 


0.27 


8 


04 


16 


41.0 


27 49 


13 


1 


.52 


± 


0.43 


43 


± 


8 


1.96 


± 


0.35 


9 


04 


16 


41.4 


27 55 


00 


1 


.31 




0.32 


22 


± 


7 


2.68 


± 


0.41 


10 


04 


16 


42.7 


27 24 


21 


1 


.02 




0.11 


58 


± 


3 


0.30 


-j- 


0.33 


11 


04 


16 


43.2 


27 40 


31 





.90 




0.11 


67 


± 


4 


1.66 




0.33 


12 


04 


16 


43.9 


28 31 


02 





.99 


± 


0.32 


—8 


± 


9 


2.75 


± 


0.31 


13 


04 


16 


50.2 


27 15 


06 


1 


.37 


± 


0.41 


35 


± 


9 


0.42 


± 


0.28 


14 


04 


16 


50.4 


27 34 


59 


2 


.66 


± 


0.50 


51 


± 


5 


1.57 


± 


0.26 


15 


04 


16 


51.3 


28 13 


14 


3 


.08 


± 


0.86 


61 


± 


8 


3.44 


± 


0.40 


16 


04 


16 


52.0 


27 36 


02 


1 


.70 


± 


0.55 


69 


± 


9 


1.28 


± 


0.28 


17 


04 


16 


53.8 


27 47 


52 


1 


.84 


± 


0.60 


45 


± 


9 


1.63 


± 


0.29 


18 


04 


16 


54.2 


27 44 


36 





.77 


± 


0.12 


36 


± 


4 


1.80 


± 


0.31 


19 


04 


16 


54.3 


27 22 


41 





.67 


± 


0.21 


51 


± 


9 


0.46 


± 


0.32 


20 


04 


16 


56.8 


27 15 


48 


1 


14 


± 


0.36 


69 


± 


9 


0.41 




0.29 


21 


04 


16 


57.2 


27 52 


46 


0.92 


± 


0.13 


25 


± 


4 


1.38 


-j- 


0.29 


22 


04 


16 


58.8 


28 08 


10 





77 


± 


0.10 


32 


± 


4 


2.23 


-j- 


0.29 


23 


04 


16 


59.7 


27 50 


29 





.87 


± 


0.14 


22 


± 


5 


1.91 


-|- 


0.26 


24 


04 


16 


59.8 


28 12 


14 


1 


.30 


± 


0.37 


49 


± 


8 


2.92 


± 


0.35 


25 


04 


17 


00.9 


28 03 


07 


1 


.07 


± 


0.27 


44 


± 


7 


2.48 


± 


0.29 


26 


04 


17 


00.9 


27 51 


43 


1 


.02 


± 


0.16 


30 


± 


4 


1.58 


± 


0.28 


27 


04 


17 


01.8 


28 21 


59 





.80 


± 


0.20 


-74 


± 


7 


5.92 




0.40 


28 


04 


17 


03.2 


28 24 


20 


2 


10 


-j- 


0.61 


— 10 


± 


8 


5.56 




0.31 


29 


04 


17 


03.8 


28 30 


26 


1 


.38 


± 


0.44 


9 


± 


9 


2.51 


± 


0.29 


30 


04 


17 


06.0 


27 48 


08 


3 


.97 


± 


0.30 


29 


± 


2 


2.11 


± 


0.28 


31 


04 


17 


06.2 


28 02 


33 


1 


.77 


± 


0.12 


41 


± 


2 


2.54 


± 


0.29 


32 


04 


17 


08.1 


27 56 


45 


1 


.24 


± 


0.14 


21 


± 


3 


2.03 


± 


0.28 


33 


04 


17 


09.7 


28 06 


09 


1 


.41 


± 


0.45 


26 


± 


9 


3.08 




0.29 


34 


04 


17 


09.8 


27 48 


51 


1 


56 


-J- 


0.48 


43 


± 


9 


2.77 




0.32 


35 


04 


17 


10.8 


28 09 


09 





.95 


± 


0.11 


22 


± 


3 


2.51 


± 


0.27 


36 


04 


17 


12.3 


27 27 46 





.99 


± 


0.15 


33 


± 


4 


1.17 




0.34 


37 


04 


17 


12.3 


27 48 


36 


3 


.51 




0.51 


34 


± 


4 


2.77 


± 


0.32 


38 


04 


17 


12.8 


27 39 


59 


2 


.32 


-|- 


0.52 


57 


± 


g 


1.81 


± 


0.32 


39 


04 


17 


13.3 


27 19 


45 





,62 


± 


0.14 


75 


± 


7 


0.52 


± 


0.32 


40 


04 


17 


16.1 


28 29 


54 


i 
i 


oo 


± 


0.52 


_3 


± 


10 


3.31 


± 


0.27 


41 


04 


17 


16.4 


28 30 


00 


1 


.34 


± 


0.20 


—8 


± 


4 


3.31 


± 


0.27 


42 


04 


17 


16.5 


28 01 


34 


1 


.41 


± 


0.10 


36 


± 


2 


2.65 


± 


0.33 


43 


04 


17 


18.7 


28 31 


11 


2 


.82 




0.47 


-6 


± 


5 


3.21 


± 


0.29 


44 


04 


17 


19.4 


27 28 


22 





.90 


-|- 


0.14 


37 


± 


5 


1.02 


± 


0.31 


45 


04 


17 


19.7 


27 18 


37 


1 


.31 


-|- 


0.41 


51 


± 


9 


1.05 


± 


0.31 


46 


04 


17 


19.9 


27 57 


09 


1 


.29 


± 


0.20 


33 


± 


4 


2.15 


± 


0.28 


47 


04 


17 


20.4 


28 29 


52 


2 


.29 


± 


0.35 


-10 


± 


4 


3.31 


± 


0.27 


48 


04 


17 


22.4 


27 57 04 


1 


.62 


± 


0.27 


31 


± 


5 


2.83 


± 


0.30 


49 


04 


17 


23.6 


27 38 


57 


1 


.76 




0.33 


70 


± 


5 


2.67 


± 


0.37 


50 


04 


17 


25.1 


27 47 


17 


2 


.19 


-|- 


0.46 


35 


± 


6 


3.70 


± 


0.41 


51 


04 


17 


25.2 


27 32 


37 


1 


11 


-|- 


0.09 


53 


± 


2 


1.57 


± 


0.29 


52 


04 


17 


25.5 


28 25 


54 


3 


.02 


± 


0.20 


1 


± 


2 


5.71 


± 


0.45 


53 


04 


17 


27.9 


28 28 


15 


3 


.23 


± 


0.20 


_7 


± 


2 


4.31 


± 


0.30 


54 


04 


17 


28.8 


27 22 


33 





.84 


± 


0.16 


47 


± 


5 


1.23 


± 


0.32 


55 


04 


17 


30.6 


28 28 


46 


3 


.25 


± 


0.44 


—5 


± 


4 


4.31 




0.30 


56 


04 


17 


31.0 


28 14 


19 


2 


.35 


± 


0.52 


40 


± 


6 


5.58 


± 


0.38 


57 


04 


17 


31.8 


27 35 


11 


1 


.86 


± 


0.40 


52 


± 


6 


3.16 


± 


0.30 


58 


04 


17 


31.9 


27 50 


13 


2 


.38 


± 


0.08 


26 


± 


1 


4.44 


± 


0.33 


59 


04 


17 


33.1 


27 11 


49 





.86 


± 


0.14 


44 


± 


5 


0.98 


± 


0.28 


60 


04 


17 


34.1 


28 30 


09 


2 


.70 




0.19 


-14 


± 


2 


4.90 


± 


0.28 


61 


04 


17 


34.5 


27 58 


49 


1 


.87 


± 


0.33 


32 


± 


5 


3.88 


± 


0.34 


62 


04 


17 


34.6 


27 44 


13 


2 


.91 


± 


0.91 


44 


± 


9 


6.67 


± 


0.39 


63 


04 


17 


34.8 


27 57 


34 


2 


.37 


± 


0.17 


28 


± 


2 


4.32 


± 


0.33 


64 


04 


17 


35.8 


27 10 


35 





.94 


± 


0.24 


53 


± 


7 


0.87 


± 


0.31 


65 


04 


17 


36.0 


28 30 


25 


2 


.67 


± 


0.62 


-15 


± 


7 


4.90 


± 


0.28 


66 


04 


17 


36.7 


27 34 


11 


2 


.43 


± 


0.31 


68 


± 


4 


2.59 


± 


0.29 


67 


04 


17 


37.5 


28 11 


23 


1 


.09 


± 


0.17 


36 


± 


4 


5.92 


± 


0.35 


68 


04 


17 


37.5 


28 09 


57 


1 


.32 


± 


0.29 


35 


± 


6 


6.32 


± 


0.36 


69 


04 


17 


37.7 


28 14 


58 


1 


.96 


± 


0.40 


21 


± 


6 


7.22 


± 


0.40 


70 


04 


17 


37.8 


27 24 


08 





.60 


± 


0.14 


36 


± 


7 


1.14 


± 


0.28 


71 


04 


17 


38.2 


28 04 


14 


2 


.10 


± 


0.28 


24 


± 


4 


4.92 


± 


0.47 


72 


04 


17 


38.6 


27 50 


47 


2 


.29 


± 


0.42 


23 


± 


5 


4.24 


± 


0.30 


73 


04 


17 


39.3 


27 15 


05 


2 


.20 


± 


0.63 


71 


± 


8 


0.95 


± 


0.25 


74 


04 


17 


39.4 


27 48 


14 


1 


.82 


± 


0.49 


35 


± 


8 


5.17 


± 


0.30 



Table 2 — Continued 



Star 




a 




8 






P 




6 






Av h 


if 


(J2000) 


(J2000) 


(%) 


(Dcg.) 




(mag 


;) 


75 


04 


17 


39.6 


28 26 


52 


4.02 


± 


0.51 


-5 


± 


4 


5.80 


± 


0.34 


76 


04 


17 


40.5 


28 09 


53 


0.97 


± 


0.15 


39 


± 


5 


6.32 


± 


0.36 


77 


04 


17 


40.9 


28 17 


16 


2.47 


± 


0.35 


7 


± 


4 


7.10 


± 


0.38 


78 


04 


17 


42.6 


28 18 


07 


2.88 


± 


0.94 


12 


± 


9 


6.49 


± 


0.49 


79 


04 


17 


43.2 


27 47 


40 


2.21 


± 


0.10 


32 


± 


1 


5.74 


± 


0.31 


80 


04 


17 


43.9 


28 05 


57 


1.72 


± 


0.23 


-4 


± 


4 


5.23 


± 


0.59 


81 


04 


17 


44.6 


27 17 


44 


1.44 


± 


0.24 


51 


± 


5 


0.93 


± 


0.26 


82 


04 


17 


45.3 


27 37 


18 


5.40 


± 


0.93 


35 


± 


5 


8.72 


± 


0.40 


83 


04 


17 


47.6 


27 32 


16 


1.81 


± 


0.52 


60 


± 


8 


5.01 


± 


0.37 


84 


04 


17 


47.6 


28 30 


30 


2.43 


± 


0.48 


-18 


± 


6 


6.26 


± 


0.29 


85 


04 


17 


47.7 


27 16 


53 


1.06 


± 


0.26 


68 


± 


7 


0.93 


± 


0.26 


86 


04 


17 


49.3 


27 52 


57 


2.61 


± 


0.53 


23 


± 


6 


4.43 


± 


0.34 


87 


04 


17 


51.1 


28 29 


16 


2.96 


± 


0.31 


-9 


± 


3 


6.72 


± 


0.31 


88 


04 


17 


51.6 


27 47 


52 


2.68 


± 


0.71 


34 


± 


8 


3.53 


± 


0.28 


89 


04 


17 


51.7 


28 26 


59 


3.81 


± 


0.38 


-12 


± 


3 


6.83 


± 


0.35 


90 


04 


17 


51.9 


28 25 


51 


3.52 


± 


0.30 


-4 


± 


2 


7.82 


± 


0.41 


91 


04 


17 


53.5 


28 26 


50 


3.42 


± 


0.53 


-6 


± 


4 


6.83 


± 


0.35 


92 


04 


17 


54.2 


27 31 


44 


2.73 


± 


0.40 


56 


± 


4 


5.01 


± 


0.37 


93 


04 


17 


55.9 


27 43 


25 


1.98 


± 


0.63 


47 


± 


9 


1.01 


± 


0.33 


94 


04 


17 


56.6 


27 27 


10 


0.93 


± 


0.22 


53 


± 


7 


3.13 


± 


0.31 


95 


04 


17 


58.8 


27 45 


00 


1.75 


± 


0.30 


49 


± 


5 


1.22 


± 


0.28 


96 


04 


18 


00.6 


27 58 


32 


1.30 


± 


0.15 


12 


± 


3 


2.41 


± 


0.30 


97 


04 


18 


00.9 


27 23 


46 


1.38 


± 


0.40 


42 


± 


8 


1.84 


± 


0.36 


98 


04 


18 


01.0 


27 21 


55 


1.13 


± 


0.32 


52 


± 


8 


1.79 


± 


0.30 


99 


04 


18 


04.6 


27 30 


01 


3.46 


± 


0.67 


57 


± 


6 


7.80 


± 


0.45 


100 


04 


18 


05.4 


28 28 


01 


3.82 


± 


0.20 


-3 


± 


1 


8.50 


± 


0.50 


101 


04 


18 


05.7 


28 22 


07 


0.67 


± 


0.13 


12 


± 


6 


12.85 


± 


0.65 


102 


04 


18 


06.3 


27 58 


51 


1.42 


± 


0.26 


13 


± 


5 


2.60 


± 


0.28 


103 


04 


18 


06.5 


27 59 


52 


0.90 


± 


0.28 


9 


± 


9 


2.71 


± 


0.27 


104 


04 


18 


06.8 


27 42 


30 


0.61 


± 


0.18 


28 


± 


8 


2.13 


± 


0.29 


105 


04 


18 


08.5 


27 50 


41 


2.13 


± 


0.21 


17 


± 


3 


1.96 


± 


0.33 


106 


04 


18 


10.5 


28 06 


43 


1.55 


± 


0.44 


3 


± 


8 


3.23 


± 


0.42 


107 


04 


18 


10.9 


28 11 


36 


1.34 


± 


0.16 


16 


± 


3 


3.61 


± 


0.38 


108 


04 


18 


11.1 


28 14 


03 


4.66 


± 


0.87 


22 


± 


5 


4.55 


± 


0.39 


109 


04 


18 


11.9 


27 53 


17 


1.51 


± 


0.46 


-2 


± 


9 


2.17 


± 


0.39 


110 


04 


18 


12.3 


27 18 


30 


0.48 


± 


0.07 


68 


± 


4 


2.25 


± 


0.37 


111 


04 


18 


13.4 


27 21 


56 


1.16 


± 


0.31 


48 


± 


8 


2.45 


± 


0.36 


112 


04 


18 


13.7 


27 41 


01 


1.19 


± 


0.15 


28 


± 


4 


2.55 


± 


0.31 


113 


04 


18 


13.9 


28 04 


56 


2.05 


± 


0.25 


16 


± 


3 


3.13 


± 


0.31 


114 


04 


18 


14.0 


28 13 


09 


1.89 


± 


0.30 


4 


± 


5 


4.55 


± 


0.39 


115 


04 


18 


15.2 


28 01 


22 


0.94 


± 


0.20 


-10 


± 


6 


2.78 


± 


0.28 


116 


04 


18 


16.8 


27 13 


02 


0.64 


± 


0.07 


33 


± 


3 


0.93 


± 


0.30 


117 


04 


18 


17.4 


27 31 


37 


3.32 


± 


0.08 


39 


± 


1 


6.97 


± 


0.40 


118 


04 


18 


22.5 


27 42 


22 


1.16 


± 


0.19 


32 


± 


5 


1.14 


± 


0.27 


119 


04 


18 


23.7 


28 08 


05 


1.02 


± 


0.22 


25 


± 


6 


2.24 


± 


0.29 


120 


04 


18 


24.8 


28 12 


27 


1.47 


± 


0.16 


19 


± 


3 


2.91 


± 


0.36 


121 


04 


18 


25.3 


27 50 


47 


1.52 


± 


0.25 


12 


± 


5 


2.29 


± 


0.37 


122 


04 


18 


26.5 


28 07 


28 


1.44 


± 


0.36 


21 


± 


7 


2.24 


± 


0.29 


123 


04 


18 


26.7 


27 16 


02 


0.78 


± 


0.17 


57 


± 


6 


1.96 


± 


0.30 


124 


04 


18 


27.3 


27 29 


07 


5.53 


± 


1.58 


36 


± 


8 


5.17 


± 


0.42 


125 


04 


18 


28.7 


27 16 


48 


1.06 


± 


0.35 


47 


± 


9 


2.79 


± 


0.28 


126 


04 


18 


30.0 


27 19 


15 


1.27 


± 


0.17 


46 


± 


4 


3.39 


± 


0.32 


127 


04 


18 


30.1 


27 54 


51 


2.09 


± 


0.40 


-2 


± 


5 


2.06 


± 


0.26 


128 


04 


18 


31.2 


27 48 


15 


1.36 


± 


0.18 


9 


± 


4 


1.58 


± 


0.36 


129 


04 


18 


34.8 


28 01 


55 


2.01 


± 


0.25 


6 


± 


4 


2.25 


± 


0.31 


130 


04 


18 


35.2 


27 11 


44 


1.38 


± 


0.44 


58 


± 


9 


2.04 


± 


0.32 


131 


04 


18 


35.3 


27 52 


21 


2.32 


± 


0.42 


18 


± 


5 


2.29 


± 


0.28 


132 


04 


18 


37.2 


27 55 


00 


2.17 


± 


0.57 


6 


± 


8 


1.91 


± 


0.26 


133 


04 


18 


38.6 


27 55 


45 


4.45 


± 


0.92 


12 


± 


6 


1.91 


± 


0.26 


134 


04 


18 


39.5 


28 02 


28 


2.09 


± 


0.19 


15 


± 


3 


2.27 


± 


0.28 


135 


04 


18 


40.1 


28 04 


28 


1.79 


± 


0.50 


4 


± 


8 


2.02 


± 


0.29 


136 


04 


18 


41.2 


27 52 


45 


1.42 


± 


0.37 


-7 


± 


7 


1.92 


± 


0.28 


137 


04 


18 


43.4 


27 33 


33 


1.14 


± 


0.30 


38 


± 


8 


1.35 


± 


0.34 


138 


04 


18 


45.2 


27 24 


49 


1.40 


± 


0.25 


19 


± 


5 


3.63 


± 


0.34 


139 


04 


18 


46.1 


27 18 


48 


3.05 


± 


0.28 


47 


± 


3 


6.93 


± 


U.oU 


140 


04 


18 


46.3 


27 21 


42 


2.85 


± 


0.34 


35 


± 


3 


7.27 


± 


0.32 


141 


04 


18 


47.3 


27 21 


29 


2.15 


± 


0.36 


29 


± 


5 


7.27 


± 


0.32 


142 


04 


18 


50.2 


27 20 


31 


2.58 


± 


0.42 


30 


± 


5 


5.49 


± 


0.31 


143 


04 


18 


53.7 


27 23 


55 


2.01 


± 


0.38 


32 


± 


5 


2.85 


± 


0.29 


144 


04 


18 


55.1 


27 20 


17 


3.54 


± 


0.47 


18 


± 


4 


7.16 


± 


0.31 


145 


04 


18 


55.9 


27 33 


00 


1.12 


± 


0.08 


25 


± 


2 


1.26 


± 


0.29 


146 


04 


18 


56.0 


27 28 


13 


0.93 


± 


0.07 


21 


± 


2 


1.70 


± 


0.29 


147 


04 


18 


57.4 


27 23 


47 


1.80 


± 


0.30 


31 


± 


5 


2.31 


± 


0.28 


148 


04 


18 


58.3 


27 48 


52 


2.27 


± 


0.60 


16 


± 


8 


1.95 


± 


0.27 


149 


04 


19 


00.2 


27 36 


09 


0.87 


± 


0.26 


29 


± 


9 


0.94 


± 


0.29 



20 



Table 2 — Continued 



Star 




,5 




P 




5 


|a 








it 


( 19000*1 












e g-) 






(mag; 


150 


04 19 00.2 


27 52 10 


2 


.02 ± 


.49 


27 


± 


7 


1 


.93 ± 0.31 


151 


04 19 05.2 


27 14 21 





.87 ± 


.11 


42 


± 


4 


13 


.04 ± 0.50 


152 


04 19 05.8 


27 21 26 


3 


.66 ± 


.88 


26 


± 


7 


2 


.14 ± 0.34 


153 


04 19 06.8 


27 21 21 


1 


.69 ± 


.26 


26 


± 


4 


2 


.14 ± 0.34 


154 


04 1 Q OK ^ 


27 17 09 


3 


.15 i 


40 


21 


-|- 


4 


7 


4- 4^ 


155 


04 19 10.0 


27 29 21 


1 


.17 ± 


.10 


20 


± 


3 


1 


.07 ± 0.29 


156 


04 19 12.1 


27 32 56 


1 


.06 ± 0. 


.11 


27 


± 


3 





.79 ± 0.25 


157 


04 19 12.2 


27 32 42 


1 


.25 ± 


.10 


16 


± 


2 





.79 ± 0.25 


158 


04 19 12.2 


27 45 31 


1 


.52 ± 


.28 


22 


± 


5 


1 


.71 ± 0.35 


159 


04 19 13.2 


27 27 33 


1 


.32 ± 


.42 


11 


± 


9 


1 


.40 ± 0.27 


160 


04 19 13.9 


27 22 59 


1 


.35 ± 


.18 


33 


± 


4 


1 


.90 ± 0.29 


161 


04 19 18.8 


27 35 24 


1 


.68 ± 0. 


.45 


36 


± 


8 


1 


.32 ± 0.23 


162 


04 19 23.1 


27 36 48 


3 


.54 ± 


.78 


24 


± 


6 


1 


.32 ± 0.23 



Note. — Right Ascension, a, is given as hours, minutes, and seconds and Declination, 8, is given as degrees, minutes, seconds. 
a Angles are equatorial, measured east from north, 
k Extinction map from Pineda et al. (2010) 



